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ABSTRACT 

We present Keck/DEIMOS spectroscopy of stars in 8 of the newly discovered ultra- faint dwarf 
galaxies around the Milky Way. We measure the velocity dispersions of Canes Venatici I, Canes 
Venatici II, Coma Berenices, Hercules, Leo IV, Leo T, Ursa Major I, and Ursa Major II from the 
velocities of 18 — 214 stars in each galaxy and find dispersions ranging from 3.3 to 7.6 km s~^. The 
6 galaxies with absolute magnitudes My < —4 are highly dark matter-dominated, with mass-to-light 
ratios approaching 1000 M0/L0,y. For the fainter galaxies, Ursa Major II and (to a lesser extent) 
Coma Berenices, we find tentative evidence for tidal disruption, which for UMa II is strongly supported 
by previous studies. If these 2 galaxies are also dark matter-dominated, they have extremely large 
mass-to-light ratios. The measured velocity dispersions of the ultra-faint dwarf galaxies are correlated 
with their luminosities, indicating that a minimum mass for luminous galactic systems may not yet 
have been reached. We also measure the metallicities of the observed stars and find that several of the 
new dwarfs have mean metallicities as low as [Fe/H] =—2.3; these galaxies represent some of the most 
metal-poor known stellar systems. The 6 brightest of the ultra-faint dwarfs extend the luminosity- 
metallicity relationship followed by more luminous dwarfs by ^ 2 orders of magnitude in luminosity. 
We detect metallicity spreads of up to 0.5 dex in several objects, suggesting multiple star formation 
epochs. UMa II and Com, despite their exceptionally low luminosities, have higher metallicities that 
suggest they may once have been much more massive. Having established the masses of the ultra-faint 
dwarfs, we re-examine the missing satellite problem. After correcting for the sky coverage of the Sloan 
Digital Sky Survey, we find that the ultra-faint dwarfs substantially alleviate the discrepancy between 
the predicted and observed numbers of satellites around the Milky Way, but there are still a factor 
of ~ 4 too few dwarf galaxies over a significant range of masses. We show that if galaxy formation in 
low-mass dark matter halos is strongly suppressed after reionization, the simulated circular velocity 
function of CDM subhalos can be brought into approximate agreement with the observed circular 
velocity function of Milky Way satellite galaxies. 

Subject headings: dark matter — galaxies: dwarf — galaxies: kinematics and dynamics — Local 
Group — techniques: radial velocities 
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1. INTRODUCTION 

The Cold Dark Matter (CDM) cosmological model 
predicts that massive galaxies such as the Milky Way 
should be surrounded by large numbers of dark matter- 
dominated satellite halos. The relatively modest pop- 
ulations of observed dwarf galaxies orbiting the Milky 
Way and Andro meda, however, seem to conflict with 
this predictio n (Kauffman. White, fc Guiderdonil 119931 : 
iKlvpin et all IT999: Moore et al. 1991 ). This apparent 
disagreement between the expected and observed num- 
bers of dwarf galaxies has become widely known as the 
"substructure" or "missing dwarf" problem. 

Proposed solutions to the substructure prob- 
lem can be broadly divided into two categories: 
cosmological and astrophysical. Cosmological so- 
lutions include mo difying the power s pe ctruni 
at small scales ([Kamionkowski fc Liddl^ 120001 : 
IZentner fc Bullocfl l2003f) and changing the proper- 
ties of the dark matter particles, such as by making 
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them warm (ICoh'n. Ayila-Reese. fc Valenzuelal 120001 : 
iBode. Ostriker. fc Turokll2001f ) or invoking a la t e deca y 
from a non-relativistic particle (jStrigari et al.l l2007aD . 
Astrophysical solutions are more prosaic, but perhaps 
easier to constrain observationally. Some of the most 
popular astrophysical solutions include the hypothesis 
that reionization could suppress the formation of dwarf 
galaxies by preventing low-mass dark matter halos from 
acquiring enough gas to form stars after z ^ 10 (e.g. , 
Bullock. Kr avtsov. fc Weinberg! 12 000: Somervillc 2005 



Benson et al. 2002 : Ricotti fc Gne din 2005; Moore et all 
20061 ) and the possibility that the dwarf galaxies we 
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observe today were once much more massive objects 
that have been reduced to the ir present a, p pearan ce 
by dramatic tidal strippiiig (iMaver et all l2001al lbl: 
iKravtsov. Gnedin. fc Klvpinll2004f) . Despite a wealth of 
ideas about how to solve the missing dwarf problem, 
distinguishing between the various proposals has proved 
to be difficult, and making sense of the tremendous 
variety of masses, luminosities, mass-to-light ratios, gas 
fractions, and star formation histories among observed 
dwarf galaxies remains a challenge. 
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Our understanding of the missing satellite problem 
and the evolution of dwarf galaxies is being rapidly re- 
vised by the discovery of a large population of new, very 
faint Local Group dwarfs in the Sloan Digital Sky Sur- 
vey (SDSS; lYork et al.ll2000h and other wide-field imag- 
ing surveys. In the past 3 years, at least 20 of these 
galaxies have been identified, nearly doubling the pre- 
viously known population. The new dwarfs inc lude 8 
additi o nal Milky Way dwarf sp heroidals ( Willma n et al.l 
2005at IZucker et al.l l2006at iB elokuro v et all l2006bl: 



Zucker et al.l l2006bl: lGrillmairl',2006; Belokurov et al I 
2007bl: ISakamoto fc Hasegawall2006f rand 1 dwarf irregu 



lar ( Irwin et al.ll2007f). 8 new d warf sph eroidals around 
Andromeda (IZucker et al]|2004[2 007: Ma rtin et al.ll2006t 
iMaiewski et al.l l2007l : llbata et al.ll2007i ). and 3 further 
new Milky Way satellites that lie in the uncertain 
paramet er space betwee n dwar f galaxies and g lobular 
clusters (IWillman et all l2005bH B elokurov et all l2007bl: 
iWalsh. Jerien. fc WillmanI 120071 ). Nearly all of these 
objects have both surface brightnesses and luminosities 
that are significantly lower than those of any previously- 
known galaxies. 

Properly placing these new discoveries within the 
framework of CDM and the missing satellite problem 
requires measurements of their internal kinematics, in 
order to determine whether the ultra-faint dwarfs are 
gravitationally bound, dark matter-dominated galax- 
ies, or tidally-disrupted systems. Only 5 of these ob- 
jects (Ursa Major I, Andromeda IX, Bootes, Canes Ve- 
natici I, and Andromeda XIV) have published stellar 
kinematics measurements, and for 2 of the 3 ultra-faint 
Milky Way dwarfs that have been s tudied already only 
a handful of st ars were o bserved (|Klevna et al.l [20051: 
Chapman et al.l [200 5: Mun oz et all l2006bl : ITbata et al l 
2006t iMaiewski et al.ii2007l ). In this paper, we present 



new stellar velocity measurements of larger samples of 
stars in 8 of the 12 new Milky Way satellites (see Ta- 
ble [1]). Including the other published and in prepara- 
tion studies that we are aware of, the only known Milky 
Way satellites that remain unobserved are Segue 1 and 
Bootes II. 

In ^J2l we describe our observations, target selection, 
and data reduction, focusing in particular on our tech- 
niques for obtaining very high precision velocity mea- 
surements with DEIMOS. We present the main results 
of this study, including measured velocity dispersions, 
masses, mass-to-light ratios, and metallicities in f|3l In 
^ we discuss the implications of our results for the Cold 
Dark Matter model and the missing satellite problem. 
We summarize our results and conclusions in ^J5l 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

We obtained spectra of individual stars in eight dwarf 
galax ies with the DEIMOS spectrograph (jFaber et al.l 
12003') on the Keck II telescope on 2007 February 12-14. 
During the observations, the weather was clear, with see- 
ing that varied between C.'S and 0'.'9 (with a very brief 
excursion to 1'.'4). The spectrograph was configured to 
cover the wavelength range 6500- 9000 A with the 1200 
£/mm grating, and the OG550 filter was used to block 
shorter wavelength light. The spectral dispersion of this 
setup is 0.33 A per pixel, and the resulting spectral reso- 



lution, taking into account our slit- width of 0'.'7 and the 
anamorphic distortion factor of 0.7, is 1.37 A FWHM 
(corresponding to 12km per pixel and 47km 
FWHM at the Ca II triplet). Exposures of Kr, Ar, Ne, 
and Xe arc lamps provided the wavelength calibration, 
and an internal quartz lamp was used for fiat-fielding. 

We observed 18 DEIMOS slit masks, with total expo- 
sure times ranging between 20 minutes and 2.5 hours. 
One to four masks were placed on each galaxy. Each 
mask contained '--^ 50 — 100 stars, of which '--^ 30 — 80% 
were expected to be actual members of the target galax- 
ies from the SDSS photometry. The positions, exposure 
times and number of slits on each mask are listed in Table 
121 Typical target stars had magnitudes of r 20 — 21. 
At r = 20, a 1 hr exposure in good seeing conditions 
yields a signal-to-noise ratio (S/N) of ~ 15, and a 2.5 hr 
exposure gives a S/N of ~ 22, where the S/N is calcu- 
lated as the average S/N per pixel in the Ca II triplet 
region. 

Target selection was carried out on star catalo gs ex- 
tract ed from the NYU-VAGC analysis (Blanto n et al.l 
I2OO5I) of th e Sloan Digital Sky Survey Da ta Release 
5 data set (jAdelman-McCarthv etlHI l2007l ).^ We set 
the target priorities so as to preferentially observe stars 
with a high likelihood of being members of the various 
dwarfs based on their color, apparent magnitude, and 
position. We constructed r, g — i color-magnitude dia- 
grams (CMDs) for each dwarf a nd overlaid globular clus- 
ter isochrones from lCleinl (|2005L hereafter C05), adjusted 
for the distance reported in the literature. We chose the 
best-fitting globular cluster red giant branch (RGB) of 
the 3 examples provided by C05. We also added a hori- 
zontal branch track derived from the M13 observations of 
C05 and asymptotic giant branch (AGB) isochrones (for 
an age of 11.2 Gyr and a metallicity of [Fe/II]= —1.3 
or Fe/H]= -1.7) from Girardi et al. (2004). The high- 
est priority targets were those located within 0.1 mag (in 
the least squares sense^) of the RGB or AGB tracks, or 
within 0.2 mag of the horizontal branch, with additional 
preference being given to brighter stars. Stars farther 
from any of the fiducial sequences were classified as lower 
priority targets. We then designed slitmasks so as to 
maximize the number of high priority targets while still 
obtaining good spatial coverage. Slitmasks were created 
using the DEIMOS dsimulator slitmask design software 
which fills in the mask area to the extent possible with 
the highest priority input targets. This automatic se- 
lection was then adjusted by hand as appropriate. The 
remaining space on the slitmasks was filled in with lower 
priority targets. The slit width for all masks was 0'.'7, and 
the minimum slit lengths were ~ 5", depending slightly 
on the density of target stars. 

In addition to the dwarf galaxy observations, we also 
obtained spectra of a radial velocity standard star, sev- 
eral telluric standards, and stars in the globular cluster 

^ The position of Leo T, which was discovered during our ob- 
serving preparations, had not yet been processed for the VAGC at 
that time, so to select targets for that galaxy we used the standard 
DR5 data. 

As in i|3.1l when we refer to the distance between a star and 
a fiducial track in a color-magnitude diagram, we mean the fol- 
lowing: dcMD = Viia - «)* - (g - i)flducial]^ -I- (r, - r-flducial)^, 
where the appropriate reference point for each star along the fidu- 
cial track is chosen so as to minimize licMD- 
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TABLE 1 

Observing Targets 



Galaxy a {J2000.0) 5 (J2000.0) My fiv'^ Distance'^ References 

(h m s) {° ' ") (mag arcsec"^) (kpc) 



Ursa Major II 


08 


51 


30.00 


63 


07 


48.0 


-3.8 


28.8 


32 


(1),{2) 


Leo T 


09 


34 


53.40 


17 


03 


05.0 


-7.1 


26.9 


417 


(3) 


Ursa Major I 


10 


34 


52.80 


51 


55 


12.0 


-5.6 


28.9 


106 


{4),(5),{7) 


Leo IV 


11 


32 


57.00 


-00 


32 


00.0 


-5.1 


28.3 


158 


(2) 


Coma Berenices 


12 


26 


59.00 


23 


54 


15.0 


-3.7 


27.4 


44 


(2) 


Canes Venatici II 


12 


57 


10.00 


34 


19 


15.0 


-4.8 


27.2 


151 


(2) 


Canes Venatici I 


13 


28 


03.50 


33 


33 


21.0 


-7.9 


28.2 


224 


(6) 


Hercules 


16 


31 


02.00 


12 


47 


30.0 


-6.0 


28.6 


138 


(2) 



References . — ( 1) fZuckcTct al. 2006b; (2) Belokurov et al. 2007H (3) llrwin et all [200^; (4) 

rWiTl man et all POOSai : (5) Belokurov et al. 2006b; (6) Zucker et al. 2006i; (7) this work 
^ Central surface brightnesses, calculated from the Plummer profile fit parameters given in the 
discovery papers cited above.'' The distances reported in the literature for these galaxies have 
generally been rounded off to the nearest multiple of 10 kpc after converting from the distance 
modulus, which is the quantity directly constrained by the data. The distances listed here have 
been calculated from the published distance moduli and rounded to the nearest kpc. 



TABLE 2 

Keck/DEIMOS Slitmask Observing Parameters 



Mask 


a (J2000) 


<5 (J2000) 


PA 


^cxp 


# of slits 


% useful 


name 


(h m s) 


(deg) 


(sec) 




spectra 


UMalLl 


08 50 38.68 


63 06 45.0 


95.8 


3600 


81 


48% 


UMan-2 


08 49 42.19 


63 11 05.6 


180.0 


3600 


87 


52% 


UMan-3 


08 53 08.75 


63 04 45.4 


109.0 


2400 


76 


62% 


UMaLl 


10 34 50.57 


51 54 47.7 


65.0 


5400 


68 


59% 


UMaL2 


10 34 22.23 


51 56 23.9 


66.0 


3600 


62 


65% 


UMaL3 


10 35 35.62 


51 56 06.4 


23.3 


5400 


68 


85% 


LeoT-1 


09 35 00.18 


17 00 56.3 


1.0 


3600 


87 


75% 


LeoIV-1 


11 32 58.69 


-00 31 41.1 


9.8 


3000 


77 


83% 


ComBer-1 


12 27 08.32 


23 52 52.0 


117.0 


9000 


78 


62% 


ComBer-2 


12 26 44.48 


23 57 58.7 


140.0 


9000 


78 


51% 


ComBer-3 


12 26 47.98 


23 54 42.8 


-20.0 


9000 


80 


65% 


CVnII-1 


12 57 12.78 


34 20 43.8 


-20.0 


9000 


67 


81% 


CVnII-2 


12 57 16.03 


34 18 51.8 


50.0 


1200 


66 


30% 


CVnLl 


13 27 59.38 


33 34 26.8 


73.0 


4140 


91 


87% 


CVnL2 


13 28 09.19 


33 31 16.0 


70.5 


4140 


94 


83% 


CVnL3 


13 28 14.34 


33 33 23.3 


-2.0 


4860 


90 


83% 


CVnL4 


13 28 02.17 


33 33 36.7 


-112.0 


9000 


115 


79% 


Herc-1 


16 31 02.70 


12 47 21.3 


104.0 


4500 


106 


83% 


Note. — Mask name, right ascension, declination, position angle and total 



exposure time for each Keck/DEIMOS slitmask. The final two columns refer to 
the total number of slitlets on each mask and the percentage of those slitlets for 
which a redshift was measured. 



NGC 1904 to serve as templates for cross-correlation with 
the dSph stars. More template observations (the globular 
cluster NGC 2419 and other radial velocity standards) 
were obtained during additional recent Keck/DEIMOS 
observing runs, with identical observing set-ups (except 
for the slit width). 

2.2. Data Reduction 

The data were reduced using version 1.1.4 of the 
DEIMOS data reduction pipeline developed for the 
DEEP2 Galaxy Redshift Survey (Cooper et al., in prepa- 
ration). Since this software was designed for faint, re- 
solved galaxies, we modified the pipeline to optimize re- 
ductions for our relatively bright unresolved stellar tar- 
gets. The main modifications were to change the cosmic 
ray rejection algorithm and to allow alignment of indi- 
vidual 2D exposures in the spatial direction before co- 
adding. In addition, we modified the long-slit pipeline to 
allow proper reduction of very bright standard stars. 



2.3. Measurement of Radial Velocities 

We measure radial velocities by cross-correlating the 
observed science spectra with a set of high signal-to- 
noise stellar templates. The stellar templates were ob- 
served with Keck/DEIMOS using the same set-up de- 
scribed above. Because template mismatch can result 
in significant velocity errors, we include a wide variety 
of stellar types and metallicities in our template library: 
giants of spectral type F8III through M8III, subgiants 
and dwarf stars. In order to cover the range of metal- 
licity expected in our low-luminosity dwarf galaxies, we 
also include several RGB and horizontal branch (HB) 
stars taken from observations of Galactic globular clus- 
ters. The stellar templates cover the metallicity range 
[Fe/H] — —2.12 to -1-0.11. All science and template spec- 
tra are rebinned onto a common wavelength array with 
logarithmic wavelength bins of size 15 km pixel"^, 
which is chosen to match the lowest spectral resolution 
present in the observed data. 
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We calculate and apply a telluric correction to each 
velocity measurement to account for velocity errors that 
result from mis-ce ntering an un resolved star within the 
slit. Following Soh n et all ()2007f ). we cross-correlate each 
science spectrum with a telluric template in the regions of 
the strong telluric absorption: 6860-6925 A, 7167-7320 A, 
7593-7690 A and 81 10-8320 A. The telluric template was 
created from the spectrum of a hot, rapidly rotating star 
(HR 1641, B3V) that was allowed to drift perpendicu- 
larly across the slit (i.e., across the 0'.'7 dimension) dur- 
ing the exposure, simulating a source that uniformly fills 
the slit, and thus accurately reflects the mean integrated 
slit function. The mean telluric offset per mask ranged 
between —7 and -1-2 km s^^, with a standard deviation 
within a mask of 3 km s^^. This correction is the veloc- 
ity error caused by the mis-centering of the science star 
within the slit from e.g., astrometry errors, small mask 
rotation, etc. Repeat observations of a number of stars 
on multiple masks demonstrate that the telluric correc- 
tion reduces the mean absolute deviation between inde- 
pendent pairs of measurements from 4.6 to 3.8 km s~^, 
reduces the weighted standard deviation of the velocity 
differences between pairs of measurements from 5.6 to 
4.2 km s~^, and improves the weighted mean difference 
from —2.0 to —0.4 km s~^, indicating that the telluric 
correction is removing both random and systematic er- 
rors from the data. 

We first calculate the telluric offset (wtou) and then 
determine radial velocities (fobs) for each science spec- 
trum. In both cases, the template and science spectra 
are continuum-subtracted; the template is then shifted 
and scaled to find the best fit in reduced-x^ space. The 
final radial velocity (v) is then: v = Vobs — vtdi — Whci, 
where Whci is the heliocentric correction determined for 
each mask. All the radial velocities presented in this 
paper include a telluric and heliocentric correction. 

The internal velocity dispersions of low-luminosity 
dwarf galaxies are, in many cases, of the same order as 
the DEIMOS velocity errors associated with individual 
measurements. In this regime, it is crucial to measure 
not only accurate velocities, but accurate velocity uncer- 
tainties. Underestimating (overestimating) the velocity 
uncertainties translates directly into larger (smaller) val- 
ues of the inferred velocity dispersion using the methods 
described in ^'d.2\ We determine our velocity error bars 
using a Monte Carlo bootstrap method, determine the 
contribution from systematic errors via repeat measure- 
ments of individual stars, and check the precision of these 
errors by comparing to higher spectral resolution data. 

For the Monte Carlo method, noise is added to each 
pixel in the one-dimensional spectrum of each science ob- 
servation based on the observed variance in that pixel. 
We assume the variance in each pixel is independent and 
distributed according to Poisson statistics. We then re- 
calculate the velocity and telluric correction for this new 
spectrum using the same routines above. Error bars are 
defined as the square root of the variance in the recov- 
ered mean velocity over 500 runs of the simulations. We 
next compare these Monte Carlo error estimates to the 
velocity differences between independent repeat measure- 
ments of individual stars. Since many of our DEIMOS 
masks covered overlapping sky areas, we placed 43 stars 
on two or more masks to obtain multiple independent 
velocity measurements. We note that one of these stars 



is likely an RR Lyrae variable (§[33]) and remove it from 
the sample of repeated observations. We are left with 
49 pairs of independent velocity measurements. The ve- 
locity difference between these independent observations 
samples the 'true' error distribution. We define a nor- 
malized error (itjv) as the velocity difference between 
two independent measurements (ui, ^2), divided by the 
quadrature sum of all error contributions: 



- /T2 — T7;2 — -272' 

where the Monte Carlo errors (ctmci, tTMC2) on each mea- 
surement are combined in quadrature with an additional 
term e equal to the error contribution from systemat- 
ics not accounted for in the Monte Carlo simulation. 
The fiN distribution should be a Gaussian of unit width. 
We therefore determine the unknown contribution from 
other systematic errors by fitting the parameter e to pro- 
duce a unit Gaussian distribution. In the left panel of 
Figure [U the final ctn distribution is plotted for the best 
fitting value of e = 2.2km s^^. The final velocity er- 
rors used in our analysis are the quadrature sum of the 
Monte Carlo and systematic errors. In the right panel 
of Figure [Ij we plot the final velocity errors as a func- 
tion of the mean per pixel signal-to-noise ratio for all 
the individual stellar velocities presented in this paper. 
Stars that fall far from the main locus of points are typ- 
ically hot horizontal branch stars (which have few sharp 
spectral features). The median velocity uncertainty for 
our sample of member stars in the ultra-faint dwarfs is 
3.4 km s~^; including mask alignment stars and bright 
foreground stars that tend to have higher S/N, the me- 
dian uncertainty for the entire data set is 2.7 km s~^. 

To demonstrate our ability to accurately measure ve- 
locities and recover velocity dispersions, we compare our 
observations of the Galactic globular cluster NGC 2419 
to higher spectral resolution Keck/HIRES observations 
of the same cluster (P. Cote 2007, private communica- 
tion). We measure radial velocities for 26 stars between 
1' and 4' of the cluster center. The HIRES data contain 
a similar number of stars in this region, although very 
few stars overlap between the two data sets. The HIRES 
spectrograph has a high spectral dispersion (0.02 A per 
pixel) and much more accurate individual velocity mea- 
surements (~ 0.95 km s""'^). We compute the reces- 
sion velocity and velocity dispersion of NGC 2419 for 
both data sets using the maximum-likelihood technique 
described in §3.21 For the recession velocity, we mea- 
sure (u)deimos = —20.7 ± 0.6 km s~^ compared to 
(u) HIRES = —21.2 ±0.5 km s~^, and for the velocity 
dispersion, ctdeimos = 2.3 ± 0.4 km s""'^ compared to 
c^HiRES = 2.3 ± 0.3km s^^. The DEIMOS observations 
agree within the 1-a limits of the more accurate HIRES 
measurements. Both sets of measure ments agree with 
the p ublished values for this cluster (jPrvor fc MevlarJ 
Il993f ). While our NGC 2419 observations have some- 
what higher signal-to-noise ratios than those typical of 
our dwarf galaxy observations, this comparison demon- 
strates that we are able to reliably measure the kine- 
matics in systems with extremely low velocity disper- 
sions (smaller than expected for the dwarf galaxies). We 
also note that we measure a metallicity based on the 
Ca triplet fines (g23) for NGC 2419 of [Fc/H]= -2.0, 
and the data are consistent with no intrinsic metallicity 
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Fig. 1. — (a) Distribution of the normalized velocity error, crjv (as defined in Eq. [T]l, for 49 pairs of repeated independent velocity 
measurements. The best fitting systematic error, e = 2.2 km s~^, is used to produce the unit Gaussian distribution shown in red. (b) 
Combined random (o"mc) and systematic velocity error for individual measurements plotted as a function of the mean per pixel signal- 
to-noise ratio. These data include all of our science targets, but not the globular cluster and standard stars observed as spectroscopic 
templates. Points that fall far from the main locus are typically hot horizontal branch stars that lack the sharp spectral features present in 
the majority of giant and dwarf stars that comprise our sample. Note that five stars in the sample have velocity uncertainties larger than 
18 km s~^, and 62 stars have S/N greater than 100, and are therefore not displayed in this plot; we chose the axis ranges so as to make 
the detailed distribution of uncertainties more visible. 



spread within the cluster. T he standard m etallicity for 
this cluster is [Fe/H]= -2.12 (lHarrislll99d) . 

As a further test of our ability to measure reliable ve- 
locities, we compare our observ ations to high r esolution 
spectroscopy in UMa I by Kle vna et al.l ()2005( ). These 
authors presented Keck/HIRES spectra of 7 stars in the 
UMa I region (5 members and 2 non-members). We re- 
observed all 7 stars with multiple measurements for 2 
stars. ^ We find excellent agreement between our mea- 
surements and those of iKlevna e t al. for 6 of the stars 
(differences within the 1 a uncertainties for 5 out of 6 
and less than 1.7 km s~^ for all 6); for star 7, both 
of our measu rements are significantly discrepant with 
Klevna et al.l's velocity. Star 7 had the lowest S/N in 
Kleyna ct aLj's observations, and they described the Ca 
triplet lines as "barely discernible above the noise." We 
conclude that either this star is a binary or variable 
star, or the velocity measured by iKlevna et al.l (120051 ) 
is in error by '--^ 8 km (their estimated uncertainty is 
5 km s-i). 

Radial velocities were successfully measured for 1015 of 
the 1460 extracted spectra across the 18 observed science 
masks. This total includes 50 duplicate measurements of 
individual stars and 124 objects identified as galaxies or 
quasars. The latter objects will be very useful as back- 
ground objects for proper motion studies and will be the 
subject of a future paper. The majority of spectra for 
which we could not measure a redshift did not have suf- 
ficient S/N. The fitted velocities are visually inspected 
to insure the reliability of the measured redshift and the 
overall quality of the spectrum. The final sample of stel- 
lar radial velocities consists of 841 unique measurements 
across the eight target dwarf galaxies. 

^ N ote that the SDSS DR,2 coordinates given in IKlevna et al.l 
II2005I ) for the target stars are up to 8" ofi' from the true positions 
as given in the DR5 data or on Palomar Sky Survey plates. 



2.4. Measurement of Equivalent Widths and 
Metallicities 

We estimate the metallicity ([Fe/H]) of individual RGB 
stars in our target galaxies using the Ca II triplet absorp- 
tion lines near A = 8500 A. We calculate the equivalent 
widths (EWs) of the three Ca II absor ption lines us- 
ing the line and continuum definitions of Rutl edge et al.l 
(jl997a). The three EWs are combined into a single 
quantity as SCa = 0.5EW(A8498A) -I- 1.0EW(A8542A) -f 
0.6EW(A8662A). We determine the error on this com- 
bined quantity with the Monte Carlo method described 
above. Added in quadrature to the Monte Carlo un- 
certainties is a systematic uncertainty of 0.3 A, which 
we determined from repeat measurements as described 
in_j^21__Wc convert SCa into m etallicity using the 
iRutledgeTl lesscr. fc Stetsoni ()1997bD empirical calibra- 
tion relationship: 

[Fe/H] = -2.66 + 0.42[I]Ca - 0.64(T/hb - V)] (2) 

The term (Vhb — V) is the magnitude difference be- 
tween the horizontal branch and the observed star, and 
corrects for surface gravity effects. We assume an ab- 
solute magnitu de for a metal-poor horizontal branch 
Mv.HB — 0.88 (|Cleml [2005). and calculate the appar- 
ent magnitude, Vhb, using the distance modulus of 
each galaxy (see Table [T|). The uncertainties in the 
distance moduli are included in the total metallicity 
uncertainties we derive. Note that assuming a single 
value for the horizontal branch magnitude in the pos- 
sible presence of multiple stellar populations may add 
an a dditional ~ 0.07 dex to the metallicity uncertain- 
ties (jKoch et al.ll2006[ ). We convert the SDSS g-band 
magnitudes into V-band using the photometric transfor- 
mations of Blanton & Rowcis (2007 ) and reddenin g cor- 
rections from | Schleg cl. Finkbcine r. fc Davis! (|1998D . The 
iRutledge et al.l calibration relation is derived for RGB 
stars in Milky Way globular clusters using the abundance 
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scale of ICarretta fc GrattonI (|1997f) . and while the cali- 
bration data only extend to [Fe/H]= —2.1 it is reasonable 
to extrapolate the relation to the slightly lower metal- 
licities found in some of our low-luminosity dwarfs (see 
i 33.4p . We restrict our metallicity analysis in i i3.4l to only 
the RGB stars in the dwarf galaxies. 

To remove foreground dwarf stars from the sample, we 
will use the equivalent width of the Na I AA8183, 8195 A 
absorption lines, which are str ongly dependent on sur- 
face gravity and tem perature (jSpinrad fc Tavlorl Il971l : 
ISchiavon et"all Il997f ). We measure the Na I equiva- 
lent width using the lin e and continuum definitions of 
ISchiavon et al.l ()1997l ). ISchiavon et al.l show that the 
Na I EW is expected to be 1 A or greater in M-type 
dwarf stars, whereas this feature i s much weaker in gi - 
ant stars at the same temperature. iGilbert et al.l (|2006l ) 
have used this feature to successfully discriminate be- 
tween dwarf and giant stars for a similar spectroscopic 
sample. 

3. RESULTS 
3.1. Selection of Members 

We use two complementary techniques to determine 
which of the observed stars are members of the dwarf 
galaxies and which are foreground stars. The first 
method classifies stars based on objective criteria: ve- 
locity, distance from the fiducial RGB and HB tracks 
(correcte d for foreground ext inction using the redden- 
ings from lSchlegel et 31111998 ' and for the distance of the 
galaxy), and the equivalent width of the Na I 8190 A 
absorption lines (described above). The exact cutoffs 
for each of these parameters needed to be adjusted in a 
few cases, but in general we use a 3 cr cutoff in veloc- 
ity (requiring a prior iteration to estimate the velocity 
dispersion), color-magnitude distance limits (defined as 
in gSH]) of 0.2 mag for RGB stars and 0.4 mag for HB 
stars, and a Na I equivalent width of less than 1.0 A. No- 
table exceptions to these cutoffs include Coma Berenices, 
which is located so nearby that we detect a number of 
subgiants, blue stragglers, and main sequence stars at 
r > 21.5 — at these faint magnitudes we extend the al- 
lowed distance from the fiducial CMD track to 0.5 mag; 
and CVn I, which has a broad giant branch that also 
necessitates widening the cutoff distance from the RGB 
track. In addition, CVn I has so many member stars 
(214) that the presence of a 3 cr outlier is likely (and in- 
deed we find one), so the velocity cutoff must be extended 
to 3.5 a, where there is only a 10% chance of finding a 
member star in our sample. 

The second method is to examine individually the fol- 
lowing properties of each star: velocity, location in the 
CMD, spatial position, fitted spectral type, metallicity, 
Na I equivalent width, and if necessary, the spectrum. 
Combining all of the available information about each 
star, and using thresholds similar to those described 
above (but less rigid), we classify each star as a likely 
member or non-member "by eye" . Both of these meth- 
ods are similar in s pirit to the techniques described by 
IGilbert et all (|2006D and iGuhathakurta et all (|2006D to 
separate M31 red giants from foreground main sequence 
stars, but without employing a full maximum-likelihood 
calculation, which is not necessary for these data be- 
cause the dSph stars are more localized in parameter 



space than M31 halo stars are, and because of the higher 
S/N. In all cases, we find excellent agreement between the 
member samples identified with the two methods (with 
occasional threshold tweaks required to produce a per- 
fect match) . A few of the galaxies contain questionable 
member stars that significantly affect the derived veloc- 
ity dispersions, and these cases will be discussed individ- 
ually in ^13. 51 We display color-magnitude diagrams and 
spatial distributions for the observed stars in each dwarf 
galaxy in Figures [2] - [D 

3.2. Central Velocity Dispersions 

Given the member samples selected in the previous 
subsection, we use the maximum-likelihood method de- 
scribed by Walker et al. (2006a) to calculate simultane- 
ously the mean velocities and velocity dispersions of each 
galaxy.'' This method assumes that the observed veloc- 
ity dispersion is the sum of the intrinsic galaxy dispersion 
and the dispersion produced by measurement errors, as 
well as that the velocity distribution is reasonably ap- 
proximated by a Gaussian. The derived velocities and 
intrinsic velocity dispersions are displayed in Table [H 
We find dispersions ranging from 3.3 ± 1.7 km s^^ for 
Leo IV to 7.6 ± 0.4 km s""'^ for Canes Venatici I. 

We plot the stellar velocity dispersions as a function 
of absolute magnitude in Figure 1 10b . There is a sig- 
nificant correlation of velocity dispersion with absolute 
magnitude, with the more luminous galaxies (My < — 6) 
having larger dispersions of ~ 7 — 8 km s~^ and the 
fainter galaxies {My ~ — 6) exhibiting smaller disper- 
sions of ~ 4 — 5 km s~^. The four low-luminosity galaxies 
Coma Berenices, CVn II, Hercules, and Leo IV are the 
first galaxies to break the velocity dispersion "barrier" 
at ~ 7 km s~^ that o bservations of the previously known 
dSphs had suggested (|Gilmore_.et al. 2007). The unprece- 
dentedly low velocity dispersions of these galaxies, and 
the correlation with absolute magnitude down to such 
low luminosities, demonstrate that if there is a floor on 
the masses of dSphs, it does not appear to have been 
reached yet. 

The likely presence of unresolved binary stars in our 
stellar velocity sample may increase the measured ve- 
locity dispers ion of our target galax ies due to binary or- 
bital motion. lOlszewski et al.l ( 19961) simulated the effect 
of binaries on the velocity dispersions of the Draco and 
Ursa Minor dSphs with very similar sample sizes and ve- 
locity uncertainties as the present study. Assuming the 
binary fractions determined for Draco and Ursa Minor 
(which range b etween .2 and 0.3 for relevant binary pe- 
riods), Olszew ski et al.l suggested that the velocity dis- 
persion from binaries alone is on the order of 1.5 km s~^. 
Since it is possible that the binary fractions may be dif- 
ferent in the lower luminosity galaxies we observed, we 
use this estimate only as a guide. For the highest velocity 
dispersion systems listed in Table [3l the effect of bina- 
ries is negligible. This result is consistent with conclu- 
sions from previous groups for other Local Group dSphs 

"* Note that the numerical values of the parameters a and b in 
Equation 9 of Walker ct al. are negative, since they are propor- 
tional to the second derivatives of Inp evaluated at the maximum 
of the function. The uncertainties on the mean velocity and the 
velocity dispersion should therefore be defined as d(u) = y^jaj and 
dcT = \/\b\ to avoid imaginary results. 
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Fig. 2. — (a) Color-magnitude diagram of observed stars in Ursa Major II. The large black circles represent stars identified as radial 
velocity members of the galaxy, the small black dots represent stars identified as non-members, and the blue crosses are spectroscopically 
confirmed background galaxies and quasars. The red curve shows the location of the red giant branch, subgiant branch, and main sequence 
turnoff populations in the globular cluster M92 and the blue curv e shows the location of the horizontal branch of M13, both corrected for 
Galactic extinction and shifted to a distance of 32 kpc (data from lClemll2005l) . (6) Spatial distribution of observed stars in Ursa Major II. 
Symbols are t he sam e as in (a) (the figure legend applies to both panels), and the ellipse represents the half-light radius of UMa II from 
IZucker et all l|2006bl V (c) Velocity histo gram of observed stars in Ursa Major II. Velocities are corrected to the heliocentric rest frame. 
The filled red histogram represents stars classified as members, and the hatched black-and-white histogram represents non-members. The 
velocity bins are 2 km s~^ wide. 
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Fig. 3. — Same as Figure[2] but for Leo T. 
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Fig. 4. — Same as Figure[2] but for Ursa Major I. 



(|Klevna et al.lll999l:IWaiker et al.ll2006aD. For our lowest correction would reduce the measured dispersion from 
dispersion system, Leo IV. the [Olszewski et al.l binary 3.3 km s^^ to 2.9 km s~^. However, this difference is sig- 
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Fig. 5. — Same as Figure[2] but for Leo IV. 
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Fig. 6. — Same as Figure [2l but for Coma Berenices. 
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Fig. 7. — Same as Figure [2l but for Canes Venatici II. 



nificantly smaller than our measurement uncertainty of 
1.7 km s~^, so we do not correct our measured disper- 
sions for the presence of binaries. Unless the binary star 
fraction in these ultra-low luminosity dwarfs is signifi- 
cantly larger than that of other dwarf galaxies, binaries 
do not significantly inflate the measured dispersions and 
inferred masses of even the lowest-dispersion dwarf galax- 
ies in our sample. 



3.3. Total Masses 

The process of determining the total mass of a dwarf 
spheroidal galaxy from the velocities of a relatively mod- 
est sample of stars that are probably located well in- 
side the virial radius of the galaxy's dark matter halo is 
fraught with difficulty. The standard technique in the lit- 
erature is to assume that (1) the galaxy is spherical, (2) 
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Fig. 8. — Same as Figure [2l but for Canes Venatici I. 
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Fig. 9. — Same as Figure [2l but for Hercules. 



TABLE 3 

Radial Velocities and Velocity Dispersions 



Galaxy 






{u)gsr 


a da Number 




(km s-i) 


(km s-i) 


(km s~^) 


(km s"-"-) (km s~-^) of stars 



Ursa Major II 


-116.5 


1.9 


-33.4 


6.7 


1.4 


20 


Leo T 


38.1 


2.0 


-58.4 


7.5 


1.6 


19 


Ursa Major I 


-55.3 


1.4 


-7.1 


7.6 


1.0 


39 


Leo IV 


132.3 


1.4 


10.1 


3.3 


1.7 


18 


Coma Berenices 


98.1 


0.9 


81.7 


4.6 


0.8 


59 


Canes Venatici II 


-128.9 


1.2 


-95.5 


4.6 


1.0 


25 


Canes Venatici I 


30.9 


0.6 


77.6 


7.6 


0.4 


214 


Hercules 


45.0 


1.1 


144.6 


5.1 


0.9 


30 



the galaxy is in dynamical equilibrium, (3) the galaxy has 
an isotropic velocity dispersion, and (4) the light distri- 
bution of the galaxy traces its mass distribution. All four 
of these assumptions may be false in reality, especially for 
the ultra-faint dwarfs that are the subject of this paper. 
SDSS photometry and foUowup imaging reveal that most 
of the dwarfs are elongated, demonstrating that they are 
not spherically symmetric systems and probably do not 
have isotropic velocity dispersion tensors. Several of the 
dwarfs also appear irregular, opening up the possibility 
that their structure has been significantly affected by the 
tidal field of the Milky Way. However, these apparently 
irregular isodensity contours could also be the result of 



the extremely low surface densities of the galaxies, which 
make their stellar distributions difficult to determine ac- 
curately. Finally, the nearly fiat velocity dispersion pro- 
files observed in all of the dSphs where spatially resolved 
kinem atics are available indicate that light does not trace 
mass ()Walker et al.ll2006aHWu II2007D . Despite these ob- 
jections, the samples of stars in the ultra-faint dwarfs 
that are spectroscopically accessible with current instru- 
ments are so small that more sophisticated analyses are 
not possible (with the exception of CVn I, which will be 
discussed in more d etail in a future pa per). We therefore 
use the method of IlllingworthI (|1976( ) to estimate total 
masses for the observed galaxies: 
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Fig. 10. — (a) Velocity dispersion as a function of absolute magnitude for the ultra-faint dwarfs. The filled black symbols represent the 
gravitationally bound dwarfs and the open gray symbol represents UMa II, which is thought to be tidally disrupted (see ^3.6| I. Circles 
are ultra-faint dwarfs in this sample and the triangle is the Bootes dSph (Martin ct al. 2007). (6) Dynamical mass as a function of total 
V-band luminosity. Symbols are the same as panel (a). The ultra-faint dwarf galaxies clearly display a trend in which the more luminous 
galaxies have larger velocity dispersions and correspondingly larger masses. Perhaps surprisingly, there appears to be a simple power-law 
relationship between mass and luminosity. 



Mtot = 167/3r,a^ (3) 

where /3 is a parameter that depends on the concentra- 
tion of the system and is generaUy assumed to be 8 for 
dSphs (.Mateo 1998), rc is the King (1962) profile core 
radius, and a is the observed central velocity dispersion. 
For most of the new dwarfs, only Plummer (half-light) 
radii rather than King core radii are available in the lit- 
erature, but we can use the fact that Tc = 0.64 * rpiummor 
to estimate the King radii. The radii and luminosities 
we have assumed for these calculations are given in Ap- 
pendix [A] Our derived masses for each galaxy are listed 
in Table [4] and plotted in Figure fTOb . We note that ob- 
jects in the bottom left corner of the plot are both the 
least massive and least luminous known galactic systems. 

The ultra-faint Milky Way satellites have masses rang- 
ing from just over 10^ (Coma Berenices) up to 
2.8 X 10'^ M0 (Canes Venatici I). Not surprisingly, CVn I, 
which is nearly as bright as previously known dSphs such 
as Ursa Minor and Draco, has a mass that is similar to 
those of the original Milky Way dSphs. Combining the 
measured masses with the absolute magnitudes listed 
in Table [U we can calculate F-band mass-to-light ra- 
tios, which are presented in Table SI The new dwarfs 
continue the trend of an anti-correlation between lumi- 
nosity aiid_^£/Tthat_has been known for many years 
(e.g.. iMateo et al]ll993^ ■ reaching mass-to-light ratios 
of ^ 1000 in l/-band solar units. Although the uncer- 
tainties on the mass-to-light ratios are substantial, owing 
primarily to the poorly known luminosities of the ultra- 
faint dwarfs, it is clear that all of these galaxies have 
quite large mass-to-light ratios. The existence of galax- 
ies with similar prope r ties to thes e was predicted recently 
bv iRicotti fc GnedinI (|2005| ) and iRead. Pontzen. fc Viell 
(|2006D . but the measured masses seem to be in better 
agreement with the models of IRicotti fc GnedinI (|2005D . 

3.4. Metallicities 

The mean stellar metallicity of a galaxy reflects the en- 
richment history of the interstellar medium at the time 
the stars were formed. We determine the mean metal- 



licity, [Fe/H], for the new dwarf galaxies based on the 
Ca II triplet equivalent width ( ij2.4p . While we can reli- 
ably measure equivalen t widths f or the majority of our 
target stars, the Rutl cdge et al.l (|1997bl ) empirical cali- 
bration that we use to convert to [Fe/H] is valid only 
for RGB stars. We therefore only include stars brighter 
than My = -fl.5 and redder than {g ~ r) > 0.3 (to 
avoid HB stars) in the metallicity analysis. We determine 
the mean metallicity and metallicity spread using the 
maximum-likelihood technique described in ^3.2\ While 
the metallicity distributions are not necessarily Gaus- 
sian, as the maximum- likelihood calculation assumes, we 
find that the mean and median of the observed metal- 
licity distributions give similar results. We run the 
maximum-likelihood algorithm twice, rejecting 3 a out- 
liers on the second run. The mean metallicities and 
metallicity spreads are listed in Table [H We find metal- 
licities ranging from [Fe/H] =-1.97 ± 0.15 for UMa II 
down to [Fe/H]=-2.31 for CVn II and Leo IV. We note 
that several of our galaxies have mean metallicities equal 
to those of the most metal-poor globular clusters and 
lower than those of other dwarf galaxies (|Harrisl [19961 : 
|Mateo"1998'), making them, along with the Bootes dSph 
(Muiioz et al. 2006b), the most metal-poor stellar sys- 
tems known. 

A strong correlation exists between the total luminos- 
ity of dwarf galaxies in the Local Group and t heir stellar 
metallicities ("Mated [19981 iGrebel et~al]l2003f l. In com- 
parison. Galacti c glob ular clusters follow no such rela- 
tionship (Harris! Il996[ ) . In Figure [11] we show that all 
but the two faintest of the ultra-low luminosity galaxies 
follow the luminosity-metallicity relationship defined by 
the more luminous dSphs. The two galaxies deviating 
from this relationship are UMa II and Com. These are 
the nearest as well as the lowest-luminosity objects in our 
sample and both galaxies (particularly UMa II) show a 
variety of evidence suggesting that they are undergoing 
tidal disruption by the Milky Way. As discussed in ? l3.6i 
we interpret the high metallicities in these two objects 
as evidence that their formation mass may have been 
significantly larger than their present mass. 

The ultra-low luminosity galaxies extend the 
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TABLE 4 

Masses, Mass-to-Light Ratios, and Metallicities 



Galaxy 


Mass 




M/Ly 


[Fc 


/H] 


0"[Fe/iI] 




(Mq 


) 




(Mq/Lq) 






Ursa Major IP 


4.9 ±2.2 


X 


10^ 


1722 ± 1226 


-1.97 


±0.15 


0.28 


Leo T 


8.2 ±3.6 


X 


10^ 


138 ± 71 


-2.29 


±0.10 


0.35 


Ursa Major I 


1.5 ±0.4 


X 


10^ 


1024 ± 636 


-2.06 


±0.10 


0.46 


Leo IV 


1.4 ± 1.5 


X 


10*5 


151 ± 177 


-2.31 


±0.10 


0.15 


Coma Berenices 


1.2 ±0.4 


X 


10^ 


448 ± 297 


-2.00 


±0.07 


0.00 


Canes Venatici II 


2.4 ± 1.1 


X 


10^ 


336 ± 240 


-2.31 


±0.12 


0.47 


Canes Venatici I 


2.7 ±0.4 


X 


10^ 


221 ± 108 


-2.09 


±0.02 


0.23 


Hercules 


7.1 ±2.6 


X 


10^ 


332 ± 221 


-2.27 


±0.07 


0.31 



UMa II may be a tidally disrupted remnant, which would artificially inflate 
its mass and mass-to-light ratio. 




Fig. 11. — Metallicity-luminosity relationship for dwarf galaxies in the Local Group. The new ultra-faint galaxies (red circles) follow the 
trend of decreasing metallicity with luminosity set by more luminous dwarf galaxies (black squares). The two lowest-luminosity objects 
(UMa II and Com) show possible evidence of tidal stripping. In comparison, Galactic globular clusters (blue triangles) do not follow any 
luminosity-metallicity relationship. Data for luminous dwarf galaxies are from Mateo (1998), Galactic globular clusters from Harris (199(|) 
and the ultra-low luminosity dwarf Bootes (open square at My = —5.8) from iMufioz et al.. (i2006b ). jMartin et al.i (2007) find a somewhat 
higher metallicity for Bootes of [Fe/H]= —2.1. The smaller horizontal bars on our galaxy measurements represent the uncertainty in the 
mean metallicity; internal metallicity spreads are indicated by the larger vertical bars. The ultra-low luminosity dwarfs are among the most 
metal-poor stellar systems in the known universe. 



luminosity-metallicity relation in the Local Group 
by an additional four magnitudes to My = —4.8. 
The location of these seven objects (including Bootes) 
on the same relationship defined by brighter dSphs 
is significant. It suggests that the stars formed in 
these galaxies are connected to the present mass of the 
galaxy and argues against significant tidal stripping, 
unless the amount of mass stripped from each galaxy 
approximately preserved the relative ordering of dwarf 
masses. We also measure significant internal metallicity 
spreads, <J[Fe/H], up to 0.5 dex in several ultra-low 
luminosity dwarfs, as listed in Table H) This suggests 
that stars formed in multiple star formation episodes, 
rather than a single burst, and firmly distinguishes 
these faint dwarfs from globular clusters, which do 
not contain mixed stellar populations. This is clearly 
the case for Leo T, which shows evidence for multiple 
stellar population s from its color-magnitude diagram 
(jlrwin et al.ll2007l ). However, a metallicity spread is the 
only evidence of multiple star formation episodes in the 
other dwarfs. Further investigation into the detailed 
abundances of these will provide a much clearer picture 



of star formation in these low mass objects. 

3.5. Comments on Individual Galaxies 

• Ursa Major II 

UMa II is one of the hardest galaxies to identify based 
on its signature in the velocity histogram (see Figure [S]!;), 
but a clear peak at —117 km s~^ emerges once the fore- 
ground dwarf stars are screened out by their Na I equiva- 
lent widths. Our measured velocity and ve locity disper- 
sion a re in good agreement with those of iMartin et al.l 
(|2007f) . We identify 20 member stars in UMa II out of 
236 targeted sources, which represents our lowest detec- 
tion rate for any of the galaxies. However, this is at least 
partly a result of our a ttempt to focus on sta rs in the out- 
lying clumps noted bv lZucker et al.l (|2006bf ) rather than 
the main body of the dwarf. There are two additional 
stars we classify as non-members that could in fact be as- 
sociated with UMa II. One of these stars has a velocity of 
—95.6 ±2.3 km s~^, just over 3 a away from the systemic 
velocity; including this star as a member would increase 
the velocity dispersion of UMa II to 8.2 ± 1.6 km s~^. 
The other candidate member, SDSSJ084947.6+630830, 
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was observed on both the first and third nights of our 
run and shows a velocity shift of 52.5 km s~^ while also 
changing spectral type from a K giant to a horizontal 
branch star. We suspect that this star is an RR Lyrae 
variable, and its mean velocity and apparent magnitude 
suggest that it is plausibly associated with UMa II, but 
we must exclude it from our velocity dispersion calcula- 
tion because of its large velocity variability. Future ob- 
servations of this star could provide improved constraints 
on the distance of UMa II. 

UMa II is a clear outlier from the My — a trend defined 
by the other galaxies in Figure 1 10b . with a dispersion of 
6.7 km s~^ despite its incredibly low luminosity. The ir- 
regular appearance, prox imity to the Milky W ay, and low 
luminosity of UMa II led lZucker et al.l (|2006b( ) to suggest 
in their discovery paper th a,t this galaxy might be i n the 
process of tidal disrup tion. iBelokurov et al.l (|2007af l and 
iFellhauer et all (|2007l ) have argued that UMa II is the 
progenitor of the recently discovered stellar tidal stream 
known as the Orphan Stream. In ? l3.6i we add addi- 
tional kinematic and abundance evidence supporting the 
hypothesis that UMa II is a tidally disrupting satellite 
that may be associated with the Orphan Stream. 

• Leo T 

Leo T is unique among the new dwarfs in that it con- 
tains gas (Mni/Afstar ~ 1) and has for med stars in the 
relatively recent past (jlrwin et al.ll2007l ). We measure a 
mean velocity of 38.1 ± 2.0 km s~^ and a stellar velocity 
dispersion of 7.5 ± 1.6 km s^^, in excellent agreement 
with the H I velocity and gas velocity dispersion mea- 
sured by Ryan- Weber et al. (in preparation). We do not 
detect any evidence for a cold stellar population to match 
the cold gas component at the center of the galaxy, but 
our sample of 19 member stars is not large enough for 
a significant detection of such a component. One would 
also only expect the youngest blue stars, which are not 
sampled by our observations, to have kinematics similar 
to the cold gas. Leo T is now one of very few dwarf 
galaxies that have well-measured kinematics from both 
the stars and the gas, and the agreement between the 
two indicates that the gas is accurately tracing the grav- 
itational potential of the galaxy. In such a small system, 
many other effects could contribute to the velocity dis- 
persion of the gas, but those contributions appear not to 
be significant. If this result also applies to other dwarfs 
then H I kinematics can be used to measure their masses 
reliably, which is useful because in many cases the gas 
extends to larger radii than the stars do. 

• Ursa Major I 

iKlevna et all (|2005[ ) reported Keck/HIRES spectra of 
5 UMa I member stars, obtaining a systemic velocity 
of —52.45 ± 4.27 km s~^ and a velocity dispersion of 
9.3li^2^ km s~^. We reobserved all 7 of the stars from 
the' Klevna et all sample (including the 2 non- members), 
and find excellent agreement on individual velocity mea- 
surements as discussed in ? I2.3I With our larger sample of 
39 member stars, our mean velocity f or UMa I is in good 
agreement with that of Kle vna et al.l . but our dispersion 
is somewhat lower than they measure (the disagreement 
is at less than 95% confidence). We rule out the ex- 
tremely high velocity dispersions of up to ~ 20 km s~^ 



allowed by iKlevna et alj's data. Usin g the same lumi- 
nosity for UMa I that IKlevna et aLl assumed (My = 
—6.75 from IWillman et al.l 1200531 ). we naturally find a 
somewhat lower mass-to-light ratio (355 ± 220 Mq/Lq) 
than they calculate, but wit h the revised magnitud e 
of My — —5.5 measured by IBelokurov et al.l (|2006bD . 
M/L becomes significantly larger. Our derived veloc- 
ity dis2ersion_Js signifi cantly lower than that obtained 
bv iMartin et all (|2007D . which may indicate that their 
uncertainties have been underestimated (see t j2.3p . We 
do not detect any evidence for the kine matically cold 
component (cr < 3.4 km s~^) suggested bv IMartin et al.l . 
despite a sample of stars that is a factor of ~ 2 larger. 

UMa I lack s a p ublished distance uncertainty. 
IWillman et al.l (|2005ar ) estimated a distance of 100 kpc 
from comparisons with the CMD of Sextans and theo- 
retical isochrones. We use fits of the M92 RGB and 
M13 HB fiducial tracks to our sample of radial velocity 
member stars to measure a more accurate distance 
modulus for UMa I of 20.13j^Q |^y mag, corresponding to 

a distance of 106lg kpc. We increase the assumed abso- 
lute magnitude of UMa I to My = —5.6 to compensate 
for this slight increase in distance. 

• Leo IV 

With only one slitmask devoted to it, and a total expo- 
sure time of less than an hour, Leo IV is the least- well 
studied galaxy in our sample. It also appears to have the 
smallest velocity dispersion, although with only 18 mem- 
ber stars and larger-than-average uncertainties on many 
of them, the uncertainty on the dispersion is significant. 
The dispersion of Leo IV also depends critically on our 
assumptions about membership. Two candidate member 
stars that we have rejected would significantly influence 
its properties if they were included. One of these stars, 
although it is located just above blue end of the horizon- 
tal branch, has a velocity oiv — 160.1±5.9 km s^^ that is 
well beyond the 3 a velocity range for the galaxy. Adding 
this star as a member would dramatically increase the ve- 
locity dispersion to 6.4 ± 2.0 km s^^. We therefore reject 
this star as a non-member (or possible binary system). 
The second star presents a more ambiguous case. It has 
a velocity of 120.2 ± 2.8 km s~^ , which is consistent with 
membership. It is located outside the half-light radius, 
but is close to 2 other member stars. However, inspec- 
tion of the spectrum reveals that the Ca triplet lines for 
this star appear to be double-peaked, and that a more 
appropriate velocity for this star may be ^ 129 km s~^. 
Given that this star may be a binary and that its true ve- 
locity is uncertain, we consider the safest approach to be 
removing it from the sample. If the star is included with 
a velocity of 120.2 km s^^ the dispersion of Leo IV would 
be 5.0 ± 1.4 km s~^; if the star's velocity is 129 km s~^ 
the effect on the velocity dispersion is negligible. 

• Coma Berenices 

Coma Berenices has the lowest luminosity of the new 
Milky Way satellites, and is located firmly in the low 
velocity dispersion half of the sample. One of the 59 
assumed member stars has a velocity just outside the 
3 cr limit, aX V — 83.1 ± 2.8 km . This star is lo- 
cated ^ 0.3 mag away from the blue edge of the sub- 
giant branch, and could be an evolved blue straggler in 
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Com. On the other hand, if we reject this star from the 
sample as a foreground contaminant, the velocity disper- 
sion of Com declines to 3.8 ±0.8 km s~^. Com is unique 
among the ultra-faint dwarfs in that its member stars 
span a wide range of g — i colors near the main sequence 
turnoff (MSTO), although only UMa II is nearby enough 
to detect such stars among the other galaxies. The pho- 
tometric uncertainties are not large enough to account 
for this spread. The stars on the blue side of the MSTO 
could be blue stragglers, but the presence of a few simi- 
larly situated stars on the red side suggests that we might 
instead be seeing the effects of multiple stellar popula- 
tions with different ages (and hence MSTO luminosities 
and colors) in Com. Although it appears in Figure ^jp 
that there may be velocity substructure in Com, this is 
partly a result of the chosen velocity binning, and a two- 
sided Kolmogorov-Smirnov (KS) test indicates that the 
observed velocity distribution is consistent with a Gaus- 
sian at the 17% confidence level. 

In terms of its luminosity, radius, and proximity to the 
Milky Way, Com is most similar to UMa II among all 
of the ultra-faint satellites. Accordingly, we consider the 
possibility that, like UMa II, Com is in the process of be- 
ing tidally disrupted. We discuss the available evidence 
for and against disruption n § 13. 6[ but conclude that Com 
is more likely to be a bound, dark matter-dominated ob- 
ject. 

• Canes Venatici II 

CVn II is a very faint and very compact dwarf with a 
low mass. We have rejected one star as a member that is 
located well outside the half-light radius and has a veloc- 
ity that is almost 4 a away from the mean velocity. If we 
instead opt to include this star, the velocity dispersion 
of the galaxy increases somewhat to 5.1 ± 1.1 km s~^. 

• Canes Venatici I 

Along with Leo T, CVn I is the only other of the new 
dwarfs to display a broad red giant branch, indicating 
a significant spread in metallicity and/or age among 
its stellar population. For Leo T, such a result is not 
surprising because young stars are visible in its CMD 
(|Irwin et al.l 120071 ) and it still contains gas, but it is 
somewhat less expected in a tiny, gas-poor system like 
CVn I. It is not currently understood how low-mass 
dSphs managed to hold on to enough gas to form stars 
over an extended peri o d of time (although see, e.g., 
iRicotti fc GnedinI 120051 iMarcohni et al.ll2006f ). but the 
same phenomenon is observed in all of the more lumi- 
nous Milky Way dSph satellites. In this sense, CVn I 
may have more in common with the previously-known 
dSphs (it approaches the lower bound of their luminosi- 
ties) than it does with its ultra-faint SDSS cousins (which 
have much lower luminosities). While other dwarfs in our 
sample show evidence for metallicity spreads, CVn I and 
Leo T — the two brightest of the ultra- faint dwarfs, and 
two of the three most massive (along with UMa I) — 
are the only two with RGB spreads as well. This result 
suggests that the critical mass and luminosity for dwarf 
galaxies to maintain the ability to form multiple stellar 
populations are My < - 7 and M > lO'^ Mq (although 
see S [3a . 

Ilbata et al.l (|2006f ) measured the velocities of 44 stars in 



CVn I and reported two distinct kinematical components 
to the galaxy: a centrally concentrated metal-rich popu- 
lation with a velocity dispersion of less than 1.9 km s~^ at 
99% significance and a more extended metal-poor com- 
ponent with a velocity dispersion of 13.912 5 C'ur 
much larger sample of 214 CVn I member stars does not 
reveal any trace of this dichotomy. Dividing our sample 
in half usin g the same metal-rich/metal-poor cutoff as 
Ilbata et al.U [Fe/H]= -2.0), we find a = 8.1±0.8 km s'^ 
for the metal-rich stars and a = 7.2 ± 0.5 km s~^ for the 
metal-poor stars. Even limiting the metal-rich selection 
to the most metal-rich ^ 10% of the sta rs (in case the 
population detected bv Ilbata et al.l[2006l represents only 
a small fraction of the overall stellar population) does 
not reveal any evidence for a cold component. We also 
do not detect any tendency for the metal-rich stars to 
be more centrally concentrated in the galaxy than the 
metal-poor stars. Finally, we measure a mean veloc- 
ity of 30.9 ± 0.6 km s^^, incon sistent with the va lue of 
- 24.5 km s-i determined bv ' Ibata et"all ()2006t ). Be- 
cause of our very large sample of member stars and our 
repeat measurements of individual stars to constrain our 
errors carefully, we conclude that there is no kinemati- 
cally cold population in CVn I, and there is no detectable 
difference between the kinematics and spatial distribu- 
tions of the metal-rich and metal-poor stars. 

• Hercules 

The only galaxy in which we detect possible evidence 
of kinematic substructure is Hercules. There is a peak 
in the velocity histogram containing 9 stars between 41 
and 43 km s^^ (compared to the mean velocity of 45.0 ± 
1.1 km s^^), with the remaining 21 stars distributed more 
broadly between 30 and 60 km s"^. Given an intrinsi- 
cally Gaussian distribution with the same mean velocity 
and velocity dispersion that we measure for Hercules, the 
likelihood of finding as many as 9 stars out of 30 in such 
a narrow, offset peak is only ~ 1%. However, a two-sided 
KS test indicates that the observed velocity distribution 
of stars in Hercules is consistent with a Gaussian at the 
43% level. We conclude that there is not yet statistically 
significant evidence of velocity substructure in Hercules. 
If the velocity substructure does turn out to be real, it 
could be a si gn in favor of th e tidal disruption hypothesis 
advanced bv lColeman et al.l (|2007D . 

Hercules contains two stars whose membership status is 
difficult to determine. The first of these, which is ex- 
cluded from our member sample, has a velocity exactly 
on the mean velocity of the galaxy, but is offset almost 
0.3 mag from the RGB and is therefore rejected on the 
basis of the photometry. Because its velocity is so close to 
the mean of the other stars, this star would have a neg- 
ligible effect on the measured velocity dispersion. The 
other candidate star, which does satisfy our membership 
criteria, has a velocity that is just within 3 cr of the mean 
velocity, at 30.7 ± 2.2 km s~^. The Na equivalent width 
and metallicity of this star are also near the edge of the 
membership ranges; if we remove this star from the sam- 
ple then the velocity dispersion of the galaxy would de- 
crease to 4.2 ± 0.9 km s^^. 

3.6. Tidal Disruption 

Two of the dwarfs presented in this paper show at least 
some evidence for on-going tidal disruption by the Milky 
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Way. As mentioned in the properties of UMa II and 
perhaps Com appear to be affected by these interactions. 

UMa II is located very close to the Milky Way, second 
only to Sagittarius (which is the archetype of tidally dis- 
rupting dwarfs) among the known dSphs. IZucker et al.l 
p006b) noted that UMa II appears irreg ular and its stars 
are bro ken up into several subclumps. iBelokurov et al.l 
()2007ar ) pointed out that the Orphan Stream lies along 
a great circle intersecting the position of UMa II, and 
our measured radial velocity of —116.5 ± 1.9km s~^ is in 
reasonable agre e ment with the 100km s^^ predicted by 
iFellhauer et all ()2007l) if UMa II is associated with the 
Orphan Stream. IFellhauer et al.l also predict a roughly 
north-south velocity gradient over several degrees within 
UMa II. Although our member sample only spans a dec- 
lination range of 13.6', we do detect a modest correlation 
between radial velocity and declination among the mem- 
ber star (correlation coefficient of —0.40), in the same 
sense as predicted. More significantly, we find strong ev- 
idence for a difference in the mean velocity between the 
eastern and western halves of the galaxy, with the stars 
on the eastern side having a velocity 8.4 ± 1.4 km s~^ 
larger than those on the western side. It is highly unlikely 
that a galaxy as small as UMa II would show significant 
coherent rotation, so this velocity gradient strongly sug- 
gests that UMa II is distorted by tidal forces. As noted 
previously, UMa II is also a clear outlier from the My — a 
trend shown in Figure 1 10b . This galaxy therefore either 
has a mass-to-light ratio several times larger than any 
other dwarf (Table H]), or its velocity dispersion has been 
inflated by the tidal field of the Milky Way. Finally, 
UMa II has a metallicity > 0.5 dex higher than would 
be expected from the luminosity-metallicity relationship 
shown in Figure 1111 Its metallicity is more appropriate 
for a system with My ~ — 10 (250 times more luminous 
than UMa II). Taken together, all of these independent 
results make a strong case for the imminent tidal disrup- 
tion of UMa II, and we are not aware of any observational 
evidence suggesting that UMa II is bound. 

Coma Berenices presents an intriguing counterpoint to 
UMa II. It shares some notable properties with UMa II, 
including an exceptionally low luminosity {My = —3.7, 
compared to My = —3.8), a location near the Milky 
Way (44 kpc instead of 32 kpc), and an unexpectedly 
high stellar metallicity. As with UMa II, we find a mod- 
est correlation of velocity with position in the galaxy 
(correlation coefficient of velocity with right ascension = 
—0.24). Dividing the galaxy in half along the minor axis, 
we find a mean velocity of 93.3±1.1 km for the north- 
western side and a mean velocity of 98.8 ±0.5 km s"^ for 
the southeastern side. This velocity difference is signifi- 
cant at the 4 cr level. As with UMa II, it is not expected 
that galaxies of this size are rotationally supported, so 
if this velocity gradient is real it suggests that Coma 
Berenices, like UMa II, may be distorted by tidal forces. 
On the other hand, there are no known tidal streams 
that are plausibly associated with Com, its velocity dis- 
persion is approximately what would be expected given 
its luminosity, and its stellar distribution is not notice- 
ably more irregular than those of the other ultra-faint 
dwarfs (there are two bright stars immediately to the 
north of Com that may be responsible for the apparent 
dis tortion of the isopleths in that direction pointed out 
bv IBelokurov et ani2007bl ). We also note that, with a 



smaller half-light radius (and larger central density; see 
M.3p than any other Local Group dwarf. Com may be 
more robust to disruption than some of its counterparts. 
While the available evidence is suggestive of the possibil- 
ity that Coma Berenices could be tidally disrupting, the 
situation is not nearly as clear-cut as it is for UMa II. We 
therefore treat Com as a bound, dark matter-dominated 
object for now, while recognizing that future observa- 
tions (most importantly, identification of an associated 
stellar stream) could change this picture. 

For the other six galaxies in our sample, we do not 
detect any statistically significant velocity gradients or 
other evidence suggesting tidal disruption. 

4. DISCUSSION 

4.1. The Missing Satellite Problem 

Understanding the nature of the ultra-faint dwarf 
galaxies and determining their impact on the missing 
satellite problem is one of the key goals of this work. 
Our observations show that with the likely exception of 
UMa II (and possibly Coma Berenices as well) the ultra- 
faint dwarfs seem to be dark matter-dominated systems, 
with masses lower than those of the previously-known 
dSphs and very large mass-to-light ratios. These galax- 
ies are currently the darkest known stellar systems in the 
universe. 

Determining the importance of the effect that the new 
dwarfs have on the abundance of satellite galaxies around 
the Milky Way requires having a way to compare ob- 
served galaxy properties to the properties of subhalos 
in N-body simulations. The simplest possible approach 
is to estimate the halo circ ular velocities of t he ultra- 
faint dwarfs as Wdrc = \/3cr (|Klvpin et al.lfl999l ). assum- 
ing that the observed dispersions are equivalent to the 
maximum dispersions reached in each galaxy, and that 
the stars have negligible orbital anisotropy. Although 
these assumptions may not be correct in detail, if we 
use cumulative satellite distributions then the results of 
this exercise are relatively insensitive to them. The cir- 
cular velocities of dark matter subhalos in the simula- 
tions can be measured robustly, giving us an appropri- 
ate point of comparison. We note that a more accurate 
means of comparing observed dwarfs to simulated sub- 
halos is to use the mass contained within 0.6 kpc, which 
is better constraine d by the observation s than the halo 
circular velocity is (|Strigari et al.|[20"07bD : these calcula- 
tions will be presented in a future paper (Strigari et al., 
in preparation) . Using the above approximation, we find 
that the ultra-faint dwarfs have circular velocities from 
i^circ = 6—13 km s^^ (for plotting and comparison pur- 
poses, we round the circular velocity of Leo IV up to 
6 km s~^). Because the fifth data release of the SDSS, 
where all of the new Milky Way satellites have been dis- 
covered, only covers 8000 deg^ of sky, we must weight 
each of the new dwarfs by a factor of ~ 5 to account for 
the additional ultra-faint dwarfs likely to be discovered 
once the rest of the sky has been similarly surveyed. 

We display the cumulative number of Milky Way satel- 
lites as a function of circular velocity in Figure[T2 We as- 
sume Poisson uncertainties on the total number of dwarfs 
(dN = VNoid + S^Nnow, where iVoid and A^ncw refer to 
the previously-known and newly-discovered dwarf galax- 
ies, respectively). For comparison, we include the sub- 
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Fig. 12. — Cumulative number of Milky Way satellite galaxies as 
a function of halo circular velocity. The filled black squares include 
the new circular velocity estimates from this paper (plus Bootes, 
but excluding UMa II) as well as all of the previously known Milky 
Way dwarfs. The open gray squares show the observed distribu- 
tion without the new ultra-faint dwarfs. We assume Poisson errors 
on the number count of satellites in each bin (computed indepen- 
dently for the new and old dwarfs), although the true uncertainties 
may be larger. The solid line plus diamonds represents the sub- 
halo abund ance within the virial radius of the Via Lactea N-body 
simulation l IDiemand et al.ll2007al) . 
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Fig. 13. — Outcome of two pro posed solutions to the missing 
satellite problem. As in Figure [T2l the filled black squares include 
the new circular velocity estimates from this paper (plus Bootes, 
but excluding UMa II) as well as all of the previously known Milky 
Way dwarfs, and the solid line plus diamonds represents the sub- 
halo abund ance within the virial radius of the Via Lactea N-body 
simulation fPiemand et al.l2007al ). The solid cyan curve shows the 
circular velocity distribution for the 51 most massive Via Lactea 
subhalos at 2 = 0. The dashed purple curve illustrates the circular 
velocity distribution for the 51 Via Lactea subhalos that had the 
largest masses at the time they were accreted by the main halo. 



halo circular velocity function from the recent Via Lactea 
simulation, currently the highest resolution (234 million 
articles) N-body simulation of a Mi lky Way-size galaxy 
Diemand. Kuhlen. fc Madaul l2007airbl ). This simulation 
assumes the best-fit WMAP three- year cosmological pa- 
rameters: flm = 0-238, »A = 0-762 , h = 0.73, n = 0.951, 
and erg = 0.74 (jSpergel et al.ll2007t ). The Via Lactea sub- 
halos include all bound halos located within the virial 
radius (389 kpc) of the main halo (see also ii4.1.2p . The 
addition of the new dwarfs, combined with the correction 
for the sky area that has yet to be observed with suffi- 
cient sensitivity, substantially changes the appearance of 
the substructure problem. The previously-known Milky 
Way satellite galaxies have a nearly flat circular velocity 
function below Ucirc = 15 km s~^, causing a discrep- 
ancy with the predictions that worsens with decreasing 
mass and reaches well over an order of magnitude below 
Wcirc = 10 km s~^. With the ultra- faint dwarfs included 
we now see a rising circular velocity function and a satel- 
lite underabundance of a factor of ~ 4 for halos with 
masses between Vdrc = 10 km s~^ and 20 km s^^. At 
fcirc = 6 km s~^ the discrepancy increases again towards 
an order of magnitude, but if the current observational 
census is still incomplete at the faint end, this is the mass 
range where that would manifest itself. The ultra-faint 
dwarfs significantly fill in the gap for satellites in the two 
lowest mass bins, but have masses that are too small to 
affect the satellite deficit at higher circular velocities. 



4.1.1. Proposed Solutions to the Missing Satellite Problem 

Using these new data, we can test a number of pro- 
posed astrophysical solutions to the missing satellite 
problem. For example, the observed dwarf galaxies could 
inhabit the most m assive subhalos at the present day 
(jStoehr et al.ll2002|). the subhalos tha t collapsed at the 
highest redshift ( Bullock et al.l |2000D , or the subhalos 



that were the most massive at the time they we re ac- 
creted by the Milky Way (jKravtsov et a l." 20041. We 
show the results of these tests in Figures [T2] and [TH To 
compare the observed dwarfs to the most massive (MM) 
subhalos, we identified the 51 halos (to match the number 
of Milky Way satellites projected to be found once the 
remainder of the sky has been surveyed) located within 
the virial radius that have the largest total masses at the 
present day in the Via Lactea simulation. The circular 
velocity function of these subhalos is plotted as the solid 
cyan curve in Figure [131 Note that because we chose the 
total number of subhalos to match the total number of 
Milky Way dwarfs, the agreement between the observed 
distribution and the cyan curve in the lowest mass bin 
is trivial. Another possibility is to compare the observed 
circular velocity function with the circular velocity func- 
tion of the subhalos that were most massive when they 
were accreted (dashed purple curve in Figure [T5)) . We se- 
lected the largest before accretion (LBA) subhalos from 
the Via Lactea simulation as the halos located within 
the virial radius of the main halo at z = that had the 
largest circular velocities at any point in the past. Again, 
the agreement at the low-mass end is simply a result of 
our choice of the top 51 subhalos from the simulation.'^ If 
the observed dwarf galaxies inhabit only the most mas- 
sive dark matter subhalos around the Milky Way, the 
shape of the mass function of the most massive subhalos 
fails to match the shape of the observed mass function. 
Using the subhalos that were most massive at the time 
they were accreted instead of the ones most massive to- 
day (i.e., allowing for mass lost by tidal stripping) brings 

^ The largest-before-accretion subsample at the present day 
(dashed purple curve in Figure I13II actually only has 46 objects 
with fcirc > 6 km s~^ because 5 of the subhalos lost so much mass 
by 2 = that they end up with even lower present-day circular 
velocities than are shown in the plot. 
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Fig. 14. — Effect of reionization on the missing satellite prob- 
lem. As in Figure 1121 the filled black squares include the new 
circular velocity estimates from this paper (plus Bootes, but ex- 
cluding UMa II) as well as all of the previously known Milky Way 
dwarfs, and the solid line plus diamonds represents the subhalo 
abundan ce within the virial rad ius of the Via Lactea N-body sim- 
ulation IIDiemand et al.l I2007al l. The solid red curve shows the 
circular velocity distribution for the 51 most massive Via Lactea 
subhalos at z = 13.6, the dashed cyan curve at 2 = 11.9, and the 
dotted blue curve at 2 = 9.6. 

the subhalo mass function slightly closer to the observed 
one, but there are still a factor of ~ 3 too few dwarfs in 
the Vcirc = 10 — 30 km s~-^ range. 

The final astrophysical solution we consider is that 
only halos that collapsed prior to reioniza tion were able 
to for m significan t numbers of stars (e.g., iBullock et al.l 
120001 : ISomervilld[200l iMoore et al.l 120061 ). Among the 
Via Lactea subhalos that are located within the virial 
radius at z = 0, we select the objects with the 51 largest 
values of Vdrc at various high redshifts. The results 
of this test are displayed in Figure [TH The solid red, 
dashed cyan, and dotted blue curves represent the sub- 
halos that would be selected if Zreion = 13.6, 11.9, and 
9.6, respectively. // reionization occurred around red- 
shift 9 — 14, and dwarf galaxy formation was strongly 
suppressed thereafter, the circular velocity function of 
Milky Way satellite galaxies approximately matches that 
of CDM subhalos. If reionization occurred at z < 8, we 
again find an underabundance of Milky Way dwarfs with 
Vcirc = 15 — 30 km s^^ compared to theoretical models, 
although we note that the individual subhalo Vcac{z) his- 
tories in the Via Lactea simulation are noisy at high red- 
shift, and the number of objects in these bins is relatively 
low. We therefore suggest that the observed mass func- 
tion of Milky Way satellite galaxies constrains reioniza- 
tion to have taken place before z = 8, in agreement with 
the 3-year WMAP results from measurement s of the cos- 
mic microwave background (zreion = 10.9"'"2;3: lPage et all 
I2006D . However, there are a number of caveats to this 
analysis: (1) the extrapolation of dwarf galaxy abun- 
dances from the SDSS DR5 sky coverage to the whole sky 
must be reasonable, (2) the observed velocity dispersions 
must provide a reasonable estimate of the halo circular 
velocities, (3) the primary physical mechanism responsi- 
ble for suppressing the formation of galaxies in low-mass 
dark matter halos must be reionization, (4) the cosmol- 
ogy used for the Via Lactea simulation (jPiemand et al.l 



l2007airb[ ) — particularly the low value of trg — must be 
a good match to the cosmology of our universe, and (5) 
the main halo simulated in Via Lactea must be a reason- 
able representation of the Milky Way. We also note that 
while WMAP and most other observational probes are 
sensitive to the mean reionization history of the universe, 
the dwarf galaxies observed in the study are sensitive 
primarily to the reionization history of the Local Group. 
If reionization was indeed responsible for the low abun- 
dance of Galactic satellites, then the Milky Way and/or 
M31 must have been undergoing vigorous enough star 
formation to ionize the intergalactic medium of the Lo- 
cal Group at z > 8. 

4.1.2. Observational Incompleteness and the Comparison 
Radius in Simulations 

One of the important assumptions involved in our anal- 
ysis in ^ M4:.l\ and 14.1.11 is the choice of the radius in the 
simulations out to which satellites should be counted. In 
the ideal case, this radius should be the virial radius, 
as we have used, but in reality the comparison between 
observations and simulations is only meaningful in the 
regime where the observations are complete. 

The observational census for Milky Way satellite galax- 
ies similar to the brighter dwarf spheroidals (My < — 9) 
should be largely complete by now; the last Milky Way 
satellite in this luminosity ra nge to be discovered was 
Sagittarius (jlbata et al.lll994 ). Recent searches of Palo- 
mar Sky Survey data, which are sensitive to such galax- 
ies anywhere within the Local Group, have not detected 
addit ional dwarfs ( Wh iting et al.l 120071 : ISimon fc Blitzl 
|2002| ) . The distribution of Milky Way dwarf galaxies as 
a function of Galactic latitude suggests that additional 
relatively bright dwarfs remain to be discovered at low 
latitudes, where extinction and foreground confusion are 
serious problems. The expected number of such objects 
if th ey are distributed uniformly arou nd the Galaxy is 
- 4 (jMateol 119981 : IWillman et"alll2004i ). which does not 
appear to be enough to significantly affect the missing 
satellite probleni; 

Very recently, iKoposov et al.l (|2007f ) have analyzed the 
detectability of faint Milky Way satellites in the SDSS 
DR5 data. They find that extremely low-luminosity ob- 
jects (My < — 5) may be missed by SDSS searches if they 
are located at relatively large distances {d > 100 kpc), as 
the horizontal branch and MSTO stars that their detec- 
tion relies on become too faint to be reliably detected 
in the SDSS. Galaxies with even lower surface bright- 
nesses than the known dwarfs (/iy > 30 mag arcsec^^), 
if they exist, are also likely to have escaped detec- 
tion. However, if there is a correlation between sur- 
face brightness and distance from the host galaxy (e.g., 
iMcConnachie fc Irwinll2006D or den sity and distance (as 
our data and iMaver et aiTl2001bl suggest), there may 
not be significant numbers of ultra-faint dwarfs at large 
distances. For more luminous dwarfs and those with 
higher central surface brightnesses, the current sample 
of Milky Way satel li tes sh ould be reasonably complete. 
The IKoposov et al.l ()2007l ) luminosity function of Milky 
Way satellite galaxies predicts that there are a total of 57 
dwarf galaxies within 280 kpc of the Milky Way over the 
whole sky, and a similar number within 420 kpc. Thus, 
our much simpler estimates of 46 ± 14 dwarf galaxies 
within 250 kpc and 51 ± 15 within 420 kpc from ^AA\ 
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appear to be well-justified. 

Incompleteness may still be a significant problem at 
the extreme faint end of the luminosity funct ion, as 
the r ecent discovery of Bootes II reveals ( Wals h~et alj 
|2007( ). Our results suggest that satellites in this lumi- 
nosity range are not gravitationally-bound dwarf galax- 
ies. If these objects are tidally disrupted dwarfs, then 
they should still contribute to the Milky Way satellite 
census, but if they are simply multiple fragments from 
larger objects (for example, if there is a physical connec- 
tion between Bootes and Bootes II) or unusual globular 
clusters then they do not correspond to dark matter sub- 
halos in the CDM simulations. Until surveys are more 
complete at faint magnitudes and some kinematic infor- 
mation is available for this class of objects, their effect 
on the missing satellite problem is not clear. 

In case there are undiscovered ultra-faint Milky Way 
dwarfs beyond d = 250 kpc, we repeated our analysis of 
using subhalos within 250 kpc from the main halo 
in the Via Lactea simulation (and discarding Leo T from 
the observed sample because it is beyond this radius). 
This smaller radius reduces the overall number of satellite 
subhalos by ^ 30%, not enough to significantly change 
our conclusions. Limiting the comparison to this radius 
shifts the preferred range of reionization redshifts slightly 
lower, but also lessens our leverage on determining the 
redshift of reionization. 

4.2. The Constant Halo Mass Hypothesis 

iMateo et al.l ()1993f ) was the first to point out that ob- 
servations of the dSphs known at that time suggested 
that they were all embedded within dark matter halos 
of mass '-^ 3 X 10^ Mq, independent of luminosity. In 
Figure [T5] we display an updated version of what has be- 
come popularly known as the "Mateo" plot, showing the 
mass-to-light ratios of all of the Local Group dSphs with 
measured kinematics as a f unction of absolute magni- 
tude. As seen previously bylMateo et al.l (|1993f ). [Mated 
(fl99M ), and iGilmore et al.l i|2007D . all of the galaxies ob- 
served prior to this work are approximate l y con sistent 
with the picture proposed by iMateo et al.l ()1993D . The 
results change, however, when the ultra-faint Milky Way 
satellites are added. Although the brightest of the ultra- 
faint dwarfs still lie within the same range of halo mass 
as their more luminous counterparts, the fainter objects 
(Hercules, Leo IV, CVn II, and Coma Berenices) are lo- 
cated well below the extrapolated trend. These galaxies 
have much lower halo masses, and hence their mass-to- 
light ratios are significantly smaller than what would be 
expected if they too were embedded in ~ 3 x 10*" Mq 
halos. Combining the new and old dwarfs, it appears 
that there are two distinct regimes: the brighter dwarfs 
(My < —9) all have similar mass dark matter halos, 
but for the fainter dwarfs {My > —9) M/L saturates 
at a value of 200 — 1000 and the halo mass declines as 
luminosity decreases (see Figure fTOb). It therefore ap- 
pears that the ultra-faint dwarfs do not occupy halos as 
massive as those of the "normal" dSphs; if there is a min- 
imum halo mass for dwarf galaxies, it is not clear that 
the observations have yet reached it. 

4.3. Central Dark Matter Densities in the Ultra-Faint 

Dwarfs 



The observed velocity dispersions and radii of the 
new dwarfs constrain their densities as well as their 
masses and mass-to-light ratios. Because these galax- 
ies are highly dark matter-dominated, the overall densi- 
ties we derive are essentially equal to the d ensity of the 
dark matter halo of each object. Following IMateo et al.l 
(|1991h . we can approximate the central density as 



(4) 



where 77 is a numerical parameter that works out to 1 
for plausible density profiles. The central densities of 
the new dwarfs range from ~ 0.08 Mq pc~^ for CVn I 
and Hercules up to ~ 2.1 Mq pc"'^ for the faintest 
galaxy. Coma Berenices. Alternatively, we can calcu- 
late mean densities from the total masses given in Table 
[4] and the core radii. To compare with the mean den- 
sities of the previously - known Milky Way dSphs tabu- 
lated by iGilmore et al.l (|2007l ). we assume that the ex- 
tent of the new dwarfs is ~ 2 King core radii. Again, 
Coma Berenices is the densest object, with a mean den- 
sity of 0.52 Mq pc-3 (= 20 GeV/c^ cm'^), almost a 
factor of five higher than any of the previously known 
dSphs. This substantially raises the limi ting mass den- 
sity o f ~ 5 GeV/c^ cm"-^ identified by IGilmore et al.l 
(|2007f ) and suggests that there may not be a true physi- 
cal ceiling on the densities of dwarf galaxies (as opposed 
to an observational ceiling) at all. However, if Com is in 
the process of tidal disruption, as our observations hint, 
then the highest dark matter densities in our sample oc- 
cur in Leo T and CVn II and are only modestly above 
the density limit of IGilmore et all (|2007l ). 

4.4. Phase-Space Density Constraints from the 
Ultra-Faint Dwarfs 

iHogan fc Dalcantonl (|2000f) introduced the parame- 
ter Q = p/a^ as an estimate of the coarse-grained 
phase-space density of t he dark matter in galaxy 
halos. As discusse d by iHogan fc Dalcantonl (2000), 
iDalcanton fc HoganI (|200H ). and IStrigari et alj (pOOQl . 
Liouville's theorem implies that observed values of Q set 
a hard lower limit on the original phase-space density of 
the dark matter. By finding the systems with the largest 
observed values of Q, we can therefore constrain the 
properties of dark matter and potentially rule out classes 
of dark matter candidates. O bserva tions of low-mass spi- 
ral galaxies by iSimon et al.l (|2005f ) yield lower lirnits of 
- 10^6 Mq pc-3 (km s-i)-3 on Q (|Strigari et al.ll2?M 
Martinez et al., in preparation), but those are less restric- 
tive constraints than are provided by the Ly-a forest. Q 
values for the ultra-faint dwarfs are listed in Table |5l 
These values are calculated under the assumption that 
the velocity dispersion of the dark matter (which is what 
Q actually depends on) is equal to the velocity disper- 
sion of the stars. Our observations show that most of the 
ultra-faint dwarfs greatly exceed the phase-space density 
constraint from the Ly-a forest, reaching a maximum of 
Q = 2.2 X 10~^ Mq pc~^ (km s~^)~^ in Coma Berenices. 
Even if the derived Q of Com has been affected by tidal 
disruption, all of the galaxies except UMa I, CVn I, and 
Hercules have Q values > 10~^ Mq pc~'^ (km s^^)^'', 
two orders of magnitude better than the Ly-a forest con- 
straint and about an order of magnitude improvement 
compared to the previously-known dSphs. In fact, the 
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Fig. 15. — Total mass-to-light ratios (in solar units) as a function of absolute magnitude for Local Group dwarf spheroidals. The red 
symbols represent the ultra-faint dwarfs from this paper (including Leo T, which is not really a dSph, and UMa II, which may be tidally dis- 
rupted, as an open red circle in the upper left). The open black squares represent all of the dSphs with previously-published kinematic data, 
including satellites of M31 as well as the Milky Way. The dashed gray lines are curves of constant dark matter halo mass (1, 2, 4, 8 X lO'^ M0 
from bottom to to p), assuming a stellar mass-to -light ratio of 2.5 Mq/Lq t/. For the previously-known Milky Way dwarfs, we recomputed 
luminosities from Irwin & Hatzidimitriou (1995) using the most up-to-date distance measurements, and then adjusted the mass-to-light ra- 
tios from the literature accordingly. References f or distance measurem ents are : Fornax (Saviane, Held, & Bertelli 2000; Mackev & Gilmor^ 

f)3; GuUieuszi k et al.ii2007 ). Leo I (BcUazzini et al.' 1 20041 ). Sculptor (Mateo'T998), Leo II (BcUazzini, Gcnnari, & Fcrraro 2005), Sextans 
30 et al. 2003), Carina IpairOra et al. 2003), Ursa Minor &ighell & Burke 1999), and Draco (Bonanos et al. 2004). References for 
M/L measurements are: Fornax (Walker ct al. 2006a), Leo I (^ohn et al. 2007; Koch et al. 2007a), Sculptor (Westfall ct al. 2006), Leo II 
jKoch et al. 2007b), Sextans (Walker et al. 2006b), Carina (Muiioz et al. 2006a), Ursa Minor (Wu 2007), Draco (Lokas, Mamon, & Prad^ 
[2005i). And II dCotiS et al.ia999i) . Cetus (.Lewis et al.ii2007.) . And IX ( Chapman et al.n2005i) . And XIV dMaiewski et al...2007i) . and Bootes 
dMufioz et al.ll2006bl) ^ 



dark matter velocity dispersion is expected to be larger 
than the stellar velocity dispersion, so the Q values we de- 
rive are upper limits on the true Q values for these galax- 
ies. Nevertheless, these Q values will further restrict the 
allowed parameter space for Warm Dark Matter parti- 
cles, and ma y have an impact on the meta-CDM sceii ario 
proposed bv lStrigari. Kaplinghat. fc Bullock! l|2007a| ). 

5. SUMMARY AND CONCLUSIONS 

We have obtained Keck/DEIMOS spectra of significant 
samples of stars in 8 of the new, ultra-faint Milky Way 
satellite galaxies recently discovered in the Sloan Digital 
Sky Survey. Using a large spectroscopic data set of radial 
velocity standard stars observed with DEIMOS, repeat 
DEIMOS measurements of stars in dwarf spheroidals and 
globular clusters, and DEIMOS and HIRES spectra of 
the same stars, we demonstrated that both our velocity 



measurements and our derived uncertainties are accu- 
rate. We then measured the velocities of 18 — 214 stars in 
each galaxy, with typical uncertainties of ~ 3.4 km . 

From our measurements of individual stellar velocities, 
we calculated velocity dispersions for each of the ultra- 
faint dwarfs. The velocity dispersions, which are listed 
in Table [3l range from 3.3 ± 1.7 km s~^ for Leo IV up 
to 7.6 ± 0.4 km s~^ for CVn I, and we showed that the 
velocity dispersions are correlated with luminosity (in- 
versely correlated with absolute magnitude). Under a set 
of simple assumptions, we calculated the total masses of 
the ultra-faint dwarfs, finding that these objects are the 
lowest-mass galaxies currently known. From the equiv- 
alent widths of the Ca triplet absorption lines we mea- 
sured the metallicities of the red giant branch stars in the 
new dwarfs and derived mean metallicities ranging from 
[Fe/H] ~ —2.0 to [Fe/H] = —2.3; several of these galaxies 
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TABLE 5 

Physical and Phase-Space Densities 



Galaxy po p Q 

(Mq pc-3) (M0 pc-3) (Mq pc-3 [km s-i]-3) 



Ursa Major IP 


1.13 


± 


0.60 


0.27 


± 


0.18 


3.7 


± 


3.1 


X 


10- 


-3 


Leo T 


0.79 


± 


0.36 


0.19 


± 


0.10 


1.9 


± 


1.5 


X 


10- 


-3 


Ursa Major I 


0.25 


± 


0.08 


0.06 


± 


0.02 


5.6 


± 


2.9 


X 


10- 


-4 


Leo IV 


0.19 


± 


0.20 


0.05 




0.05 


5.3 


± 


9.9 


X 


10- 


-3 


Coma Berenices 


2.09 


± 


0.86 


0.52 


± 


0.24 


2.2 


± 


1.4 


X 


10- 


-2 


Canes Venatici II 


0.49 


± 


0.25 


0.12 


± 


0.07 


5.1 


± 


4.1 


X 


10- 


-3 


Canes Venatici I 


0.08 


± 


0.02 


0.02 


± 


0.01 


1.7 


± 


0.5 


X 


10- 


-4 


Hercules 


0.10 


± 


0.04 


0.02 


± 


0.01 


7.7 


± 


5.2 


X 


10- 


-4 



UMa II may be a tidally disrupted remnant, which would artificially 
inflate its density. 



are the most metal-poor stellar systems yet discovered. 

We summarize our primary conclusions from this study 
as follows: 

1. The ultra- faint Milky Way satellites are dark 
matter-dominated dwarf galaxies with lower 
masses than any other known galaxies. 

2. The only clear exception among the 8 galaxies we 
observed, as well as those previously observed by 
others, is Ursa Major II Base d on its clumpy ap- 
pearance (jZucker et alj l2006bl), s mall galactocen- 
tric distance fZucker et al."2006b) , associated tidal 
stream (Fcllhauer ct al. 2007, inflated ve- 
locity dispersion fi j3.2p . unusually high metallic- 
ity ( ij3.4p . and possible velocity gradient ft j3.5p . we 
conclude that UMa II is in the late stages of tidal 
disruption. The other dwarf with My ~ — 4, 
Coma Berenices, has a similarly high metallic- 
ity that indicates it may have suffered substantial 
tidal stripping as well. Because Com lacks most 
of the other supporting evidence for tidal disrup- 
tion, we assume for now that it is still a bound, 
dark matter-dominated dwarf, although we recog- 
nize that future observations may show otherwise. 
Based on these results, we suggest that My ~ —4 
(3.4 X 10^ L0) is the lower limit to the luminosity 
of gravitationally-bound dwarf galaxies. We there- 
fore predict that objects such as Willman 1, Segue 
1, and Bootes II will prove to be tidally-disrupted 
remnants. 

3. The 6 ultra- faint dwarfs with My ~ — 4 follow 
the luminosity-metallicity relationship established 
by the more luminous Local Group dwarfs, and ex- 
tend the relation by ~ 2 orders of magnitude in 
luminosity. The faintest dwarfs, UMa II and Com, 
are both outliers from this relationship, with metal- 
licities more than 0.5 dex too large for their lumi- 
nosities (or conversely, luminosities that are more 
than two orders of magnitude too small for their 
metallicities). We detect metallicity spreads of up 
to 0.5 dex in several objects, suggesting multiple 
star formation epochs. 

4. The total mass-to-light ratios of the ultra-faint 
dwarfs reach as high as 1000 Mo/LQ^y (UMa I). 
While the brighter galaxies {My ~ — 9) have mass- 
to-light ratios consistent with the hypothesis that 



all dwarf spheroidals are embedded within dark 
matter halos of the same mass, the fainter galaxies 
depart from this trend in the sense that their mass- 
to-light ratios are too low (i.e., they have lower 
masses). We therefore suggest that the minimum 
mass for dwarf galaxies (as opposed to the mini- 
mum luminosity mentioned earlier), if there is one, 
may not have been reached yet. 

5. The ultra- faint Milky Way satellites, after correct- 
ing for the sky area not covered by DR5 of the 
Sloan survey, substantially increase the abundance 
of dwarf galaxies with very low masses (wcuc < 
15 km s~^), thereby reducing the satellite deficit 
compared to CDM simulations to a factor of ^ 4. 
Proposals to remedy the missing satellite problem 
by placing the observed dwarf galaxies in the most 
massive dark matter subhalos (at the present day) 
around the Milky Way or in the subhalos that were 
most massive at the time they were accreted by the 
Milky Way do not reproduce the observed shape 
of the circular velocity function. If we assume in- 
stead that only the halos that acquired a signifi- 
cant amount of mass (ucirc ~ 8 km s~^, varying 
somewhat with Zreion) before the redshift of reion- 
ization were able to form stars, t hen the subha- 
los fro m the Via Lactea simulation (jPiemand et al.l 
l2007ai rbl) approximately match both the total num- 
ber of Milky Way dwarfs and the shape of the cir- 
cular velocity function. 

6. The central dark matter densities of the ultra-faint 
dwarfs are as high as 2.1 Mq pc~^ (0.8 Mq pc~^ 
if Coma Berenices is tidally disrupting), signifi- 
cantly larger than those of the previously-known 
dwarf spheroidals. The phase-space densities are 
also higher than those of other astrophysical sys- 
tems (Q > 10"^ Mq pc-3 (km s"i)"3), which wiU 
place significant limits on non-CDM dark matter 
models. 



Data presented herein were obtained at the W. M. 
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nership among the California Institute of Technology, 
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tics and Space Administration. The Observatory was 
made possible by the generous financial support of the 
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TABLE 6 

Parameters of the Ultra-Faint Milky Way Satellites 



Galaxy 


Mv 




distance 


Plummer 


Plummer 


References 








(kpc) 


(arcmin) 


(pc) 




Ursa Major II 


—3.8 ± 


0.6 


32I4 


13.6 


127 ± 21 


{1):(2) 


Leo T 


-7.1 ± 


0.3 


417t?o 


1.4 


170 ± 15 


(3) 


Ursa Major I 


-5.6 ± 


0.6 


106+g 


10.0 


308 ± 32 


(2), (4), (5) 


Leo IV 


-5.1 ± 


0.6 


158tii 


3.3 


152 ± 17 


(6) 


Coma Berenices 


-3.7 ± 


0.6 


44 ±4 


5.0 


64 ± 7 


(6) 


Canes Venatici II 


-4.8 ± 


0.6 


i5it;^ 


3.0 


132 ± 16 


(6) 


Canes Venatici I 


-7.9 ± 


0.5 




8.5 ±0.5 


554 ± 63 


(7) 


Hercules 


-6.0 ± 


0.6 


138^^3 


8.0 


321 ± 36 


(6) 


Segue 1 


-3.0 ± 


0.6 


23 ± 2 


4.5 


30 ±3 


(6) 


Willman 1 


-2.5 ± 


1.0 


38 ± 7 


1.9 ±0.3 


21 ± 5 


(8) 


Bootes II 


-3.1 ± 


1.1 


60 ±6 


4.1 ± 1.6 


72 ±28 


(9) 


Bootes 


-5.8 ± 


0.5 


60 ±6 


13.0 ±0.7 


227 ± 26 


(3) 



References. — (1) IZucker et al."2006bl: (2) D. Zucker & V. Belokurov 2007, 
private c ommunication; (3 ) Irwin et al. 2007; (4) Belokurov ct al. 2006t j; (5) this 
wo rk; (e) IBelokurov eraLii2007h ; (7') , Zucker et al.ii2006a.; (81 .Willman et al.||200l ; 
(Ql IWalsh et al.ll2007l 
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APPENDIX 

APPENDIX: COLLECTING UNIFORM DATA FOR THE ULTRA-FAINT DWARF GALAXIES 

The new ultra- faint dwarf galaxies have been discovered by a number of different authors, which means that their 
properties have not all been determined in a uniform manner. In order to calculate masses and mass-to-light ratios 
as consistently as possible for each galaxy, and to facilitate future studies of these objects, in this appendix we collect 
the currently available data on all 12 of the new Milky Way satellites (see Table [5]) ■ 

With the exception of Willman 1 and Bootes II, we use the absolute magnitudes and radii determined by the 
Cambridge group (Zucker, Belokurov, Irwin et al.) from SDS S data. This includes the revised absolute magnitude 
of My — —5.5 for UMa I reported in iBelokurov et "all (|2006bf ). which differs substantially from the original value of 
My — —6.75 estimated bv lWillman et al.l ()2005aD . although the uncertainties on both numbers are admittedly large. 
Note that our improved distance for UMa I of 106 kpc (compared to the pre viously r eported 100 kpc) requires a 
corresponding change in the absolute magnitude to My = —5.6. The Belokuro v et al.l magnitude for UMa I is not 
accompa nied by an uncert ainty; by analogy to the other galaxies of similar luminosity we assume an uncertainty of 
0.6 mag. IBelokurov et al.l have also re-measured the radius of UMa I in the same manner as they did for the other 
Milky Way galaxies, finding a Plummer radius of 10' (D. Zucker & V. Belokurov 2007, private communication). 

UMa II is described in the discovery paper only as having an angular extent of 0?5 x 0?25 and a half-light radius 
of approximately 120 pc (Zucker et al. 2006b). Similarly to UMa I, the authors have re-fit the light profile using the 
same method as they did for the other new dwarfs and measured a half-light radius of 13.'6 (D. Zucker & V. Belokurov 
2007, private communication). Note that unlike the Plummer radii we use for the other dwarfs, this radius is the mean 
of the Plummer radius and the exponential scale radius; in most cases the two radii are very similar. 

Bootes II does not have a published distance uncertainty, so given the angular proxi mity to Bootes and the ap parently 
identical distance moduli, we assume the same distance uncertainty for Bootes II as IBelokurov et all ()2006bl ) derived 
for Bootes. 

Most of the new discoveries do not have published uncertainties for their Plummer radius fits. For these objects, we 
assume an uncertainty of 6% (the uncertainty given for CVn I) on the angular radius for the purposes of calculating 
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the uncertainty on the corresponding physical radius. 
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